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Durante la compactacion hidraulicamente inducida de una 
capa granular surgen patrones de fractura. Con simulaciones 
numericas estudiamos como estos patrones dependen de las 
propiedades del gas intersticial asi como de las propiedades 
del medio poroso. En particular aquf se estudia en detalle la 
relacion entre la velocidad de propagation de fractura y la 
presion de inyeccion. 



During the hydraulically induced compaction of a granular 
layer fracture patterns arise. In numerical simulations we 
study how these patterns depend on the gas properties as 
well as on the properties of the porous medium. In particular 
the relation between the speed of fracture propagation and 
injection pressure is here studied in detail. 
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INTRODUCTION 

Stress induced by fluid or gases can cause diverse materi- 
als to break and fracture. Such hydraulic fractures are a 
natural and common phenomenon in the field of volcan- 
ism and are artificially initiated to enhance the recovery 
of natural gas and mineral oil by fracturing the reservoir 
rock with pressurized fluids. Recently a new perspective 
on hydrofractures was added with the storage of supercrit- 
ical CO2 attracting the interest of an increasing number 
of researchers. In this respect two scenarios are consid- 
ered. First it is one option to inject CO2 into existing 
hydrofractures, and second the injection of the CO2 can 
create additional fractures [1, 2]. The typical components 
for such fractures are a porous material and a compressible 
gas. Injection of pressurized fluids in a porous material, 
deforming beyond the elastic limit, has been studied in 
granular materials in Hele-Shaw cells, [3, 4, 5, 6, 7, 8], with 
the injection of air or oil in systems with open boundary 
conditions, and during cyclic loading [9]. It was also stud- 
ied in systems with a confinement for the grains, prevented 
from getting out of the cell, which allowed to observe the 
formation of thin fractures [10]. In this paper [10] it was 
found and discussed a criterion that the porous media and 
the fluid need to fulfill to allow the formation of fractures. 
For this purpose the gas' viscosity was varied. It was fur- 
ther discussed how the shape of the fractures depend on 
the properties of the porous material and of the injected 
gas in simple 2 dimensional (2D) numerical simulations. 
In contrast to the previous article we will not change the 
properties of the injected gas or the porous material in this 
present article. Here we explore in particular the effect of 
the amplitude of the fluid pressure imposed in the source 
on the fracture morphology. Furthermore, all simulations 



here will be ran in a regime where fractures are created. 

SIMULATION SETUP 

As shown in Fig. 1 the setup consists of a cell with two 
glass plates separated by 1 mm. The gap between the 
plates is filled with particles. The empty space between 
the grains is saturated with a fluid that has the same prop- 
erties as the fluid that is injected. Consequently, the only 
two media involved in the dynamics are the grains and 
the fluid. At the start of the simulations the average solid 



volume fraction of the grains is p 



(o) 



0.42. This start- 



ing solid volume fraction is homogeneous with negligible 
density fluctuations although the particles are at random 
positions. The value of pi°^ = 0.42 is chosen to be less 
than the possible maximum ofp[ max) =0.60 to allow com- 
paction of the grains. On the inlet side of the cell the 
pressure is imposed to Pj. Several simulations are per- 
formed with Pj ranging from a value of Pj = 0.5 • 10 5 Pa 
to a value of Pj = 2.5 • 10 5 Pa above the atmospheric pres- 
sure of Po = 1-0 • 10 5 Pa. On the opposing side to the inlet 
the cell has an open boundary for the fluid but particles 
are not able to leave the cell here. In a real experiment, 
this could be achieved by using a net with a mesh smaller 
than the particles. The remaining boundaries are com- 
pletely sealed for both media. In the simulation around 
200 000 grains of diameter 140±10% pm are involved. Fi- 
nally, the pressure at the inlet is increased and maintained 
as a step function in time, at a steep ramp, and particles 
hardly move before the maximum injection pressure Pj is 
reached. 
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Figura 1. Numerical setup of the system. 

THEORY AND MODEL 

Using a well tested numerical model we have the freedom 
to explore the parameter space independently. The details 
of the method can be found in [6, 11, 12, 13, 14, 15, 16, 17], 
and alternative models can be found in [18, 19, 20, 21, 22]. 
The model describes the fluid in terms of a pressure field 
while the porous medium is modeled by simulating discrete 
particles. 

THEORY AND SIMULATIONS 

Dynamics of the gas phase: The equation for the evolu- 
tion of the pressure P = P + Po, where Po is the atmo- 
spheric pressure and P the local pressure fluctuations is 
given by: 



dP 
~dt 
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In this equation the pressure is described in terms of the 
local granular velocity u, the viscosity pf of the gas, the 
local porosity (j) = 1 — p s and the local permeability n. 
Eq. 1 is derived from mass conservation of the gas and 
the granular medium and by assuming a local Darcy law. 



Dynamics of the particles: 
use Newton's second law: 



For the particles we basically 



m—r- = Fi 

dt 
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plates. 



is the contact area of the particles with the 



RESULTS 

We ran a set of six simulations for injection pressures of 
Pi = (0.5, 1.0, 1.5, 2.0, 2.5, 3, 0) • 10 5 Pa above atmospheric 
pressure Po, a fluid viscosity of pf = 18.0 mPa-s and a 
friction coefficient with the side plates of 7A = 4.0. The 
injected gas is considered as an ideal gas and has the com- 
pressibility of air /3t = 1 /Po at Po . The value of Pj at the 
inlet is reached very fast and particles start to move shortly 
after. During this compression of the particles fractures 
emerge in the granular packing. In Fig. 2 a set of snap- 
shots of the particle density is shown. The snapshots are 
taken at increasing time from left to right. Each horizontal 
row of pictures corresponds to one of the six simulations at 
a different injection pressures. In Fig. 3 snapshots of the 
corresponding pressure field in the cell are displayed. The 
pressure field is normalized to one to allow a qualitative 
comparison. 
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with particle velocity v p , particle mass m, particle mass 
density p m , cell spacing h and the number density p n = 
p s pm/m. Fi are linear inter-particle solid contact forces. 
Fd is a viscous damping force during particle collisions. 
For P a , the interaction with the side plates we assume 
that the normal stress P^ in the granular packing is pro- 
portional to the in-plane stress p]' by a factor A (Janssen 
hypothesis). Using further a Coulomb friction model we 
state that the frictional force F a per particle with the glass 
plates is proportional to the normal stress by a friction 
coefficient 7. With these two assumptions we find an ex- 
pression for the friction force with the side plates. 



F a < 7 S a (2P q ± ^p m gh) 



■jS a (2\p\\ 



■ Pmgh). 
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Figura 2. Snapshots during the simulations of the particle density 
in the Hele-Shaw cell, displayed for decreasing injection pressure Pj 
from top to bottom and as a function of time (left to right). Low 
particle density appears brighter in the snapshots. Under air 
injection, fractures and fingers of low particle density emerge and 
propagate, x- and y-axis units are given in cm. The y-axis specifies 
the distance from the inlet. 
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Figura 3. (Color online) The pressure evolution for decreasing 
injection pressure Pj (top to bottom) and as a function of time (left 
to right). High pressure appears yellow (brighter) in the snapshots. 
x- and y-axis units are given in cm. 

In these plots an apparent feature is the different propaga- 
tion speed and position of the emerging fractures. A high 
injection pressure causes the fractures to propagate faster. 
To quantify this observation we can plot the position of 
the most advanced finger tip as a function in time. This 
is done in Fig. 4. The plot clearly proves the previous 
observation. Furthermore it turns out that the systematic 
increase of the propagation speed is also proportional to 
the square root of the injection pressure. This is checked 
in Fig. 5. Here the rescaling of the fracture tip position 
by the square root of the injection pressure \[Pl results 
in a collapse of the graphs. The disagreement at the later 
stages of the simulations in this plot results from the finite 
size of the system, which allows fractures to grow only up 
to a certain size. Finally we can state that the fingers grow 
according to: 



Y t = VP~if(t). 



(4) 
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Figura 4. (Color online) The position of the most advanced 
finger/fracture as a function of time at different injection pressure 
Pj. The higher the injection pressure Pj the further fingers grow. 




Where f(t) is a function which appears in the plots to 
be almost linear at early stages of the finger growth for 
t < 0.01s. 
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Figura 5. (Color online) The position of the most advanced 
finger/fracture rescaled by the square root of the injection pressure 
\[P[. As a function of time the graphs at different injection 
pressure Pj collapse onto a single graph. 

In Fig. 2 we also observe that the fingers at high injection 
pressure propagate further into the packing before com- 
plete compaction of the grains takes place. This can be 
also seen in Fig. 4 where the finger position stops growing 
at longer distances from the inlet the higher the injection 
speed is. 

Apart from the finger position, the increase of the injec- 
tion pressure also affects the shape of the fingers. In Fig. 
2 it can be seen that the fingers get more branched and 
fracture-like at higher injection pressure. At low injec- 
tion pressure fingers appear to be straighter while increas- 
ing the injection pressure, fingers develop more and more 
branches. At the highest injection pressure of Pj = 3, 0-10 5 
Pa the fingers clearly show characteristics of fractures. 
At high injection pressure Pj the pressure gradients are 
the largest. When the boundary is deformed the expected 
changes of the pressure gradients are therefore also higher 
at high injection pressure than at low injection pressure. 
At low P/, one expects a lower pressure gradient every- 
where, and thus a low effect of seepage forces and a slower 
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deformation. Leading to overall smoother pressure gradi- 
ents and a more stable front deformation. We thus expect 
faster finger propagation, and more branching at a higher 
injection pressure. 

Finally also the spatial distance between the fingers de- 
pends on the injection pressure (see 6). At low injection 
pressure the number of fingers is higher than at high injec- 
tion pressure, as can be also seen in Fig. 2. In general the 
finger spatial frequency decreases in time after injection 
has started and fingers propagate through the cell. This 
can be shown by calculating the average of the character- 
istic spatial finger wavenumber in x direction. First the 
powerspectrum Sj of each horizontal line j of the particle 
density is calculated. Taking the average of these power 
spectra results in a single power spectrum S. From this av- 
erage distribution of wave numbers the characteristic wave 
number (k) in x-direction is defined and calculated in the 
following way, using an average of k with the power spectra 
as a weight: 



(k) 
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(5) 



The results in Fig. 6 show a decrease of finger frequency 
in time. As a trend we notice that at higher injection 
pressure the finger frequency decreases faster than at low 
injection pressure. However the simulation at P/=250 
kPa differs from the other simulations. In this simulation 
we also observe a finger propagating directly along the 
right boundary in Fig. 2. Close to the wall this finger 
appears to propagate faster than the other fingers in this 
simulation. Because the simulation at Pj =250 kPa is the 
only simulation where this appears it also stands out in 
the plots for the average wave number Fig. 6. This is 
presumably due to a finite size effect, and such outlier is 
frequently met in granular systems, which are known to 
present a large variability and sensitivity on details of the 
initial state, (see e.g. [17]). Otherwise for higher injec- 
tion pressure does the finger frequency not only decreases 
faster but also drops to a lower value before the grains get 
compacted. 

This coarsening of the finger frequency is the result of 
two mechanisms. First the pressure gradient between the 
finger tip and the outlet gets higher the closer the finger 
tip moves to the outlet. Assuming a linear pressure profile 
though the porous media the pressure would drop to zero 
on a shorter and shorter distance the closer the finger 
advances to the outlet. At the same time the gas also 
leaks into the side walls of the finger. This increases the 
pressure in the porous material around a finger. In the ar- 
eas where this pressure increase takes place less advanced 
neighboring fingers would thus experience a lower pressure 
gradient. The speed of this fingers is thus reduced. This 
means the more a finger advances to the outlet the faster it 
moves. At the same time the pressure increase in the area 
around an advanced finger decreases the pressure gradient 
in front of less advanced fingers. This causes the less ad- 
vanced fingers to propagate slower or to stop completely. 
This mechanism will result in a coarsening of the finger 
frequency. Further more we expect this mechanism to be 
active on a typical length scale which is comparable to the 



skin depth of the pressure profile. This mechanism also 
appears in the basic Saffman Taylor instability [23] . 
A second mechanism that will account for a coarsening of 
the finger frequency is the compaction of the grains on the 
sides of a finger. During the propagation the finger width 
increases and branches at a 90 degree angle arise on the 
sides of fingers. This compacts the granular material on 
the sides of an advancing finger. How far this compaction 
propagates on the sides depends on the properties of the 
granular material and also on the finger width and how the 
side branches develop. Where this compaction has occur 
preceding fingers are slowed down or stopped. The size 
of the compaction front around the fingers sets a second 
length scale for the coarsening of the finger frequency. 
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Figura 6. (Color online) The average of the spatial finger 
wavenumber in x direction. 

CONCLUSIONS 

The increase of the injection pressure primarily causes fin- 
gers to propagate faster through the granular packing. 
Fingers at high injection pressure also tend to be more 
branched and fracture-like than the fingers at low injec- 
tion pressure. It was shown that the position of the frac- 
ture propagation in time increases with the square root 
of the injection pressure \fPj. Furthermore we discussed 
the observed coarsening of the characteristic spatial fin- 
ger wave numbers in terms of two mechanisms. A first 
mechanism that controls the coarsening arises from the 
fluid seepage into the granular media. Where the length 
scaled for this mechanism was argued to be of the size 
of the pressure skin depth. To further explain the coars- 
ening a second mechanism causing the coarsening of the 
finger wavenumber was highlighted. This second mech- 
anism introduces a length scale for the coarsening with 
the size of the compaction front in the granular mate- 
rial around a finger. Acknowledgments: We thank Gus- 
tavo Sanchez-Colina for help about the Spanish grammar. 
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